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Abstract 

We provide an explicit rigorous derivation of a diffusion limit - a stochas- 
tic differential equation with additive noise - from a deterministic skew- 
product flow. This flow is assumed to exhibit time-scale separation and has 
the form of a slowly evolving system driven by a fast chaotic flow. Under 
mild assumptions on the fast flow, we prove convergence to a stochastic dif- 
ferential equation as the time-scale separation grows. In contrast to existing 
work, we do not require the flow to have good mixing properties. As a con- 
sequence, our results incorporate a large class of fast flows, including the 
classical Lorenz equations. 



1 Introduction 

There is considerable interest in understanding how stochastic behaviour emerges 
from deterministic systems, both in the mathematics and applications literature. 
In this note we provide a simple explicit construction of such emergent stochastic 
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behaviour in the setting of skew-product flows exhibiting time-scale separation. We 
prove a diffusion limit for the following ordinary differential equations (ODEs): 

± (e) =e-7o(y (e) ) + /(x (e) ,y (e) ), x (e) (0)=£, (1.1a) 
y^=e- 2 g(y^), y®(0) = 77. (1.1b) 

Here G X, y^ G ^ with A" = M rf and ^ = Roughly speaking we assume 

that the equation for y( £ ) has a compact attractor A C y supporting an invariant 
measure and satisfying certain "mild chaoticity" assumptions. These conditions are 
stated precisely in Assumptions 11.21 below. In addition, we assume that /o should 
average to zero with respect to this measure. 

Consider the stochastic differential equation (SDE) 

X{t) = £ + / F{X{s)) ds + VxW{t). (1.2) 
J 

where W is unit <i-dimensional Brownian motion, E is a d x d covariance matrix 
(depending on /o and g) and F(x) is the average of f(x, •) with respect to the 
aforementioned invariant measure. The goal of the note is to prove the following 

limit theorem relating the solution of (11.11 ) and X of (11.2I ). Throughout we use 
—} w to denote weak convergence in the sense of probability measures fflUl- 

Theorem 1.1. Let Assumptions \L2\ hold and let r/ be a random variable distributed 
according to the measure \x on y and fix any £ G X. Then, almost surely with 
respect to 7] and W, there is a unique solution (x, y) G C 1 ([0, 00); X x y) of (11.11 ) 
and a unique solution X G C([0, 00); X) of (11.2I ). Furthermore x^' -^- w X in 
C '([0, 00), X) as e^O. 

Throughout the note we make the following standing assumptions. 

Assumptions 1.2. The differential equations (11.11) satisfy the following: 

1. The vector fields g : y — >• and fo'.y^-W 1 are locally Lipschitz, and 
the vector field f : X x y — >■ R d is bounded and Lipschitz with uniform 
Lipschitz constant L. 

2. Equation (ll.lbl ) vviY/z e = 1 /las a compact invariant set A, r/ G A, a«<i 
f/jere jj arc invariant probability measure fx supported on A; expectation with 
respect to this measure is denoted by E. 

3. The vector field fo'.y—t averages to under fi : E/o = 0. 

4. Define W n {t) = /„"* / (y (1) (r)) dr, for t > 0. Fix any T > 0. Wfe 

assume the weak invariance principle (WIP), namely that W n —)- w V^W 
in C([0, T], M. d ) as n — > 00 for unit d-dimensional Brownian motion W and 
some covariance matrix E, independent ofT. 
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5. Define 



F(x) = Ef(x, •) = f f(x,y)n(dy) 



J A 



noting that it is globally bounded by |/|oo an d globally Lipschitz. with con- 
stant L, since fi is a probability measure. Let T > 0. We assume the follow- 
ing large deviation principle (LDP): 



where b(a, T) is independent ofx and b(a, T) ->0asT-> oo for all a > 0. 

Remark 1.3. (a) The regularity conditions on f,fo,g in assum ption 1 guara ntee 
global existence and uniqueness of solutions to the ODEs (11.11 ) and the SDE (11.21 ) 
for all positive time and all initial conditions £ 6 X, rj £ A. We note that unifor- 
mity of the Lipschitz constant for f is automatic in y since A is compact and we 
are at liberty to modify f outside A. 

(b) Assumptions 4 and 5 hold for a large class of flows. In particular, the WIP 
and LDP are proved in ft3\ |?]/ for flows that have a Poincare map modelled by a 
Young tower ^5] [6]/ with summable tails. This includes Anosov and Axiom A flows, 
nonuniformly hyperbolic flows such as Henon-like flows ( where the Poincare map 
has a Henon-like attractor), and Lorenz attractors ^ (including the case of the 
classical parameter values in [8]). In this class of examples, the Poincare map has 
good statistical properties and limit laws such as the WIP and LDP transfer to the 

flow m. 

There are two main routes leading to emergent stochastic behaviour in deter- 
ministic systems. The first is through the elimination of a large number of degrees 
of freedom, and the reliance on the central limit theorem to provide fluctuations 
and the second is through time-scale separation; see flOl for an overview. The 
first mechanism does not require any assumption of chaotic behaviour and may be 
observed in large systems of oscillators; work in this area was initiated in ifTTl and 
more recent work includes |[T2l [131 . The second mechanism relies on the pres- 
ence of some fast chaotic dynamics to induce white noise and has a lon g history 
in the applied literature; we mention, in particular, the work in lfl4l [151 1161 1171 . 
Our work provides a very simple scenario in which the second mechanism may be 
used, provably, to establish emergent stochastic dynamics. We anticipate that the 
basic ideas would apply to a far larger class of problems as indicated, for example, 
by the program outlined in [18]. Moreover the basic mechanism that underlies the 
work in this note was identified and studied in the seminal paper [19]]. However 
the conditions in that paper can be hard to verify for specific ordinary differential 
equations. In contrast our construction holds for explicit systems on y such as the 
classical Lorenz equations. 

An important aspect of our theory is that we require no knowledge of mixing 
properties of the flow, In contrast, previous rigorous results in the literature required 
strong assumptions on the mixing properties of the flow. See [20] for the most 
powerful results in this direction where it is required that the flow has stretched 
exponential decay of correlations. Even for Anosov flows this has been proved 
only in very special cases |[2Tll22ll23l. Su perpolvnomial decay has been proved for 
typical Anosov and Axiom A flows [24 , 25 ] and typical nonuniformly hyperbolic 
flows governed by Young towers 11261 12711 but only for very smooth observables; 
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this smoothness would have to be imposed on /o. For the Lorenz equations there 
are currently no results at all on rates of mixing (though superpolynomial decay 
holds for typical nearby flows by j27l ). 



2 Weak Invariance and Large Deviation Principles 

The WIP and LDP assumed above have two useful implications which we detail in 
this section, and then use in the next section to prove the main theorem. 

Proposition 2.1. Let [x^ e > (t) , y^ e \t)) denote the solution to (11.11) with f = 0, 
£ = and with r\ a random variable distributed according to the measure [i on 
A. Let T > 0. Then x^ -+ w V^W in C([0,T],X) as e -> 0. Here, W is unit 
d-dimensional Brownian motion and the covariance matrix £ is independent ofT. 

Proof. Note that y^\t) is the solution to the IVP y = g(y), y(0) = r\. Define 
W n (t) = n~l J n * My {1 \r)) dr, for t G [0,T]. By the WIP, W n -> w V^W in 
C([0,T],X) asra oo. 

Now y^(t) = y^(te~ 2 ). Hence 

xW(t) = e- 1 f f (y®(s))d8 = e t f (yW '(r)) dr. 



Writing n = e 2 , we obtain x^ e \t) = W n (t) and the result follows. □ 
Proposition 2.2. For any a,T > 0, n > 0, x G M d , 
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■ f{x,yW(r))dT - F(x) < a + 2\f\ 00 b(a,T). 

JnT 

Proof. Let S(T, x) = |^ f Q T /(x, y (1) (r)) fir - F(x)\. Then, by the LDP, 
ES(T, x) = [ S(T, x)dfi+ I S{T, x) dfi 

JS(T,x)<a JS(T,x)>a 

< a + |S|oo 6(a,T) < a + 2\f\ oa b(a,T). 
This proves the result for n = 0, and the general case follows from invariance of 



3 Diffusion Limit 

We now prove the diffusion limit contained in Theorem ll.il The method of proof 
generalizes that described in Chapter 18 of |[28l for homogenization in SDEs with 
additive noise and a skew-product form. 
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Proof of Theorem 17771 To prove weak convergence on [0, oo), it suffices to establish 
weak convergence on [0, T] for each fixed T > 0. 

Write W^(t) = j* * 7/0(2/^ ds. By integrating the x^ equation we have 

x (e \t)=t+ [ lf (yU( 8 ))da+ f f(xV(s),yV(s))ds 
Jo € Jo 

= i + W^ ) {t) + Z i - e \t)+ [ F(x {e) (s))ds 

Jo 

where 

Z (e) (t) = f {f(x& (s) , (*)) - F(x® (*))) cfa. 
Jo 

We show below that Z& -> Oin L 1 (C([0, T], AT); //). By Proposition O -►, 
in C([0,T], AT). It follows that W< e ) + -+ w \fcW in C([0,r], X). 
Now consider the continuous map Q : C([0, T],A^) — > C([0,T],X) given by 
G(u) = v where v is the unique solution to the integral equation 

v(t) =£ + u{t)+ f F(v(s))ds. 
Jo 

Then 

x ( e ) = g(w {t) + z®). 

Since continuous maps preserve weak convergence, it follows that x^ — X as 
required. 

It remains to show the L 1 convergence of Z^ to in C([0,T],X). De- 
fine g(x,y) = f(x,y) - F(x) and note that < 2|/|oo and Lvp(g) < 2L. 
Then Z^(t) = f* g(x^(s),y^(s)) ds. Let N = [t/6] and write Z^(t) = 
Z&(N6) + Io where I„ = j t NS g(x^(s),y (e \s)) ds. We have 

|/o| < (t-NSMvo <2\f\ 00 S. (3.1) 

We now estimate Z^\N$) as follows: 

W-l /•(n+l)<5 

z( £ >(JV5) = £ / ^ (e) (s),y (e) 00)cfa 

n=0 • /n<5 
W-l /.(n+l)5 

= E / (5(x (£) («),y (£) (s)) - 5 (x (£) (n5),yW( S ))) 

n=0 • /n5 

W-l /-(n+l)<5 

+ E / 5(x (£) (^),y {e) (^))^ 

n=0 ' /n < 5 

= /l + / 2 - 
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For s e [nS, (n + 1)5], we have \x^(s) - x^\n5)\ < (|/ |oo + |/|oo)*e _1 . 
Hence 

\h\ < iV<5Lip( 5 )(|/o|oo + l/loo)^" 1 < 2L(|/ | oo + \fUT8e~ 1 . (3.2) 
Next, 

N-l f (n+l)S N-l 
n=0 JnS n=0 

where 

J n = J- 1 / f{x^{n5),y^(s))ds- F(x&(n5)) 

J nS 

f (n+l)Se~ 2 

= 8~ l e 2 I f(x®(n8),yW(s))ds-F(x®(n5)). 

J nSe~ 2 

Let a > 0. Applying Proposition 12.21 E|J n | < a + 2|/|oo 6(a, 5e~ 2 ). It follows 
that 

E sup \I 2 {t)\<T{a + 2\f\ 00 b{a,5e- 2 )). (3.3) 

*6[0,T] 

3 

Finally, we set 5 = e?. By (13.11) and (13.21) . lim^o E sup tg [ 0T ] |io(*)l = 
lim e _ + oIEsup tg [o i r] I -MO I = 0- Moreover, by (13.31 ) and the assumption on b in 
the LDP in Assumption 1 1.21 5). we have 

limsupE sup 1^2(^)1 < Ta, 
e->0 te[o,T] 

and, since a > is arbitrary, 

limE sup \I 2 (t)\ = 0. 
te[o,T] 

Altogether, lim e _s.o E sup 4g [ jT i |Z( e )(0| = as required. 



4 Conclusions 

The construction in this paper shows how some new ideas in the theory of dynam- 
ical systems can be used to prove a homogenization principle in ODEs, leading to 
emergent stochastic behaviour. The arguments are very straightforward, and are 
given only in the case of additive noise. However in the situation where the lim- 
iting SDE is one dimensional the ideas of Sussman |29l can be used to derive a 
limiting SDE in which noise appears multiplicatively. Generalizing these ideas to 
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skew product flows where the SDE is of higher dimension will require the theory 
of rough paths [30 ] and is the subject of ongoing work. 

Finally a comment on the diffe rences between j28l homo geni zation and av- 
eraging in ODE systems like dl.ll) . There is current interest [31] in the deriva- 
tion of averaging principles for systems of ODEs exhibiting three time scales of 
order 0(e -2 ), C(e _1 ) and 0(1). The motivation is the construction of efficient 
numerical sc heme s for computation of the averaged solution, which is determinis- 
tic. Theorem ll.il which also concerns the limiting behaviour of a system contain- 
ing three time-scales, corresponds to a homogenization principle with a stochastic 
limit, rather than an averaging principle with determinitic limit. Thus our work 
provides an example of a three scale system for which an effective deterministic 
averaged equation cannot exist. 
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